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Abstract 

Taylor-Goldstein equation ( TGE ) governs the stability of a shear-flow of an 
inviscid fluid of variable density . Using a canonical class of its transformations 
it is investigated from a geometrical point of view . 

Rayleigh's point of inflection criterion and Fj0rtoft's condition of instability 
of a homogenous shear-flow have been generalized here so that only the profile 
carrying the point of inflection is modified by the variation of density . This 
fulfils a persistent expectation in the literature . 

A pair of bounds exists such that in any unstable flow the flow-curvature ( a 
function of flow-layers ) exceeds the upper bound at some flow-layer and falls 
below the lower bound at a higher layer . This is the main result proved here . 

Bounds are obtained on the growth rate and the wave numbers of unstable 
modes , in fulfillment of longstanding predictions of Howard . A result of Drazin 
and Howard on the boundedness of the wave numbers is generalized to TGE . 

The results above hold if the local Richardson number does not exceed 1 / 4 , 
otherwise a weakening of the conditions necessary for instability is seen . 

Conditions for the propagation of neutrally stable waves and bounds on the 
phase speeds of destabihzing waves are obtained . It is also shown that the set 
of complex wave velocities of normal modes of an arbitrary flow is bounded . 

Fundamental solutions of TGE are obtained and their smoothness is 
examined. Finally sufficient conditions for instability are suggested. 
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1 introduction 



The object of our study here is Taylor-Goldstein equation (TGE). It was dis- 
covered by G. I.Taylor and S.Goldstein , independently and simultaneously , in 
connection with their studies on ' parallel shear-flow' of a fluid of variable density 
( see Taylor 1931; Goldstein 1931 ) . It governs the stability of the flow . TGE is 
important because of its applications to oceanography and meteorology and 
has received considerable attention in the literature on hydrodynamic stability 
( see Drazin & Howard 1966; Drazin & Reid 1981 ) . It describes waves in the 
ocean and clouds in the sky . Here we explore the fascinating mathematical 
structure . that this equation is naturally associated with . 

In the limit of vanishing density variation TGE collapses into the classical 
equation of Rayleigh ( 1880 ) which has a well developed theory ( see Lin 1955 ) 
characterized by its elegance , close agreement with experiments and geomet- 
rical flavour . The fundamental result of this theory is the celebrated point of 
inflection criterion of Rayleigh for the instability of a parallel shear-flow . 

Theory of shear-flows has many beautiful results but much remains to be 
understood yet . The hterature on the subject ( see Drazin & Howard 1966; 
Drazin & Reid 1981 ) is full of expectations that under suitable conditions the 
characteristic features of Rayleigh-theory will carry over to TGE . There 
also are laments to the effect that the point of inflection loses its significance 
in the context of TGE . In spite of persistent and valuable attempts starting 
with Synge ( 1933 ) , Yih ( 1957 ) , Drazin ( 1958 ) up to S . Friedlander ( 2001 ) 
( see Drazin & Reid 1981 ; Friedlander 2001 ) , these expectations have found 
limited fulfillment . There also are standing questions . The predictions of 
Howard (1961) are yet to be settled and it is not known if the well known 
condition (Richardson number is less than 1/4 somewhere in the flow) of 
Howard ( 1961 ) and Miles ( 1961 ) is sufficient to ensure instabihty ( see Fried- 
lander 2001 ) . Further there is no simple method for solving TGE . Thus there 
is a need for a detailed investigation . 

To understand TGE , we construct a large class of its transformations which 
transform it into a canonical form resembling Rayleigh's equation and develop a 
method that yields information on stability via each transformation of the class . 
The entire information is then focused upon the question of stability and made 
independent of transformations and coherent . This results in a description of 
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the inner core of TGE . A central role in this description is seen to be played 
by a characteristic function of the flow , we call the flow-curvature . 

This analysis establishes the point of inflection criterion for TGE, general- 
izes results from Rayleigh-theory , shows that the condition of Howard and 
Miles does not ensure instability , proves the predictions of Howard ( 1961 ) , 
relates instability to the crossing of a pair of bounds in a deflnite order by the 
flow-curvature , examines the propagation of neutral and unstable waves and 
yields fundamental solutions of TGE . 

Taylor-Goldstein equation is w"{z) + A{z) w{z) = (1-1) 
with boundary conditions w{zi) = = w{z2) . (1-2) 

Here zi < z < Z2 , w{z) is a complex-valued C^— function on 
[^1,^2] 5 the prime " ' " denotes differentiation d/dz and 



where g > is a constant, /3{z) = — p'{z)/po > (unless 
stated otherwise ) , the flow velocity u{z) and the fluid density p{z) are 
sufficiently smooth real-valued functions , the average density po is a positive 
constant and c = Cr + ici is a complex constant called the complex wave 
velocity . A; > is the wave number , q the phase speed and k Ci the 
growth rate of a harmonic perturbation wave ( A; , c ) . 

The stability of ' parallel shear-flows ' is an important problem in classical 
theory of hydrodynamics . Many eminent investigators have contributed to its 
understanding ( see Friedlander & Howard 1998 ) . TGE governs the stability 
of a 2-dimensional parallel shear-flow of an inhomogeneous fluid . 

We imagine an inviscid and incompressible fluid of variable density flowing 
in the (x, z)-plane, between the lines z = zi and z = Z2 . Here the 
axis is vertical , directed upwards and the x— axis is horizontal i being 
the unit vector along it . 

Suppose at the point ( x , z ) the velocity of the fluid particle is u{z) i 
and the density is p{z) . Thus a layer of the fluid characterized by a fixed 
value of z has a constant velocity and density . It experiences a force due to 
buoyancy, proportional to the density gradient P{z) , exerted on it by the 
neighbouring layers . These incompressible and parallel layers of the fluid slip 



A{z) = -k^ - 



u"{z) 
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smoothly on each other and generate a shearing movement in the fluid . This 
time-independent flow is cafled a paraflel shear-flow and denoted by {u ,/3) . 

The ( hnear ) instabihty of a shear-flow is perceived of in the theory of sta- 
biUty , as the existence of a harmonic perturbation wave , which vanishes at the 
boundaries , whose amphtude grows exponentiafly with time and whose prop- 
agation is permitted by the hnearized form of the equation of motion of the 
flow — a nonhnear PDE due to Eiiler . 

Only stable parallel shear-flows are seen to persist in the laboratory or in 
nature . Governed by the equation of motion , the unstable flows even if they 
be parallel to start with , soon become time-dependent and nonparallel . 
This is because small and time-dependent natural perturbations grow with time 
in these flows . Thus stability is directly related to the question : which parallel 
flows persist ? 

Here we adapt from classical hydrodynamics the well known criterion for the 
stability of this flow (see Drazin & Reid 1981 ), as a deflnition as follows : 

A harmonic perturbation wave {k , c) is called a normal mode of the 
flow {u,/3) if equations (1) and (2) have a nontrivial smooth solution w{z) . 
A normal mode is called unstable if q > and neutral if q = . A 
flow is called unstable if it has at least one unstable normal mode and stable 
otherwise. 

The most fundamental result on TGE was conjectured by Taylor in 
1931 and proved by Howard ( 1961) and Miles ( 1961 ) . Howard transformed 
TGE to prove that for any unstable flow {u, (3) , the flow discriminant A{z) 
is positive for some z in {zi, Z2) , where A{z) = [ u''^{z) — 4g/3{z) ] . 
Equivalently if Richardson number [ g /3{z) / u''^{z) ] > 1/4 everywhere, 
then the flow is stable . 

Synge (1933) proved that Wmin < < Wmax for any unstable mode 
{k , c) of a flow {u, /3) . Howard ( 1961 ) went deeper and elegantly proved 
that in this case the wave velocity c must lie in the upper half plane , inside 
the semicircle whose diameter is the real interval [wmin, ^max] • Kochar &; 
Jain ( 1979 ) took a step further and replaced the semicircle by a semi-ellipse 
inside it with the same base . Howard also proved that the growth rate k Ci 
of unstable modes is bounded above by V ^max / 2 , generalizing a result due 
to Holland ( 1953 ) ( see Drazin & Howard 1966 ) . Further he predicted that 
kci ^ as A;— ^00 ( this is proved easily for homogeneous flows ) and 
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also that under suitable conditions the wave numbers k of unstable modes of 
a flow are bounded . 

The special case of TGE when u = and (3 takes arbitrary real values 
was known to Rayleigh . In this case instability is expected due to gravitational 
overturning if a heavier layer of the fluid lies above a lighter layer . Rayleigh 
provided the mathematical support for this observation and proved that the 
static flow {0 , (3) is unstable if and only if /3{z) < for some z . 

We now discuss the basic results on homogeneous flows {u , ) . When 
(3 = 0, (1.1) is called Rayleigh's equation. Rayleigh (1880) proved the 
fundamental result that in any unstable flow, u"{z) must assume both 
positive and negative values so that the flow- velocity u{z) must have a 
point of inflection at some point Zg in {zi, Z2) . Fj0rtoft ( 1950 ) added that 
in this case u"{z){u{z) — u{zs)} < for some z. 

Friedrichs ( 1942 ) ( see Drazin & Howard 1966 ) using Sturm- Liouville 
theory proved the existence of a smooth , neutrally stable mode ( ks , u{Zs) ) 
if for some Zs in [2:1,2:2], u"{zs) = 0, K{z) = — u"{z) / {u{z) — u{zs)} 
is integrable and K{z) > tt^ 7(2:2 — 2:1)^ whenever u{z) u{zs) . Here 



the minimum being taken over all functions such that and 0' are 
square integrable and (f){zi) = = ^(2:2) . 

Drazin & Howard ( 1966 ) proved that if K{z) is integrable and K{z) > 
whenever u{z) u{zs), then k'^ < k^ for any unstable mode {k, c) . Their 
argument also shows that K{z) > [tt^ 7(2:2 — 2:1)^ + k"^] somewhere. The 
last statement follows also from independent arguments given by M.B.Banerjee 
et.al. (2000). Thus if < K{z) < 'k'^/{z2-ziY whenever u{z) ^ u{zs) , 
then the flow is stable . 

Heisenberg ( 1924 ) and Tollmien ( 1929 ) ( see Drazin & Reid 1981 ) obtained 
neutral wave solutions ( Cj = ) of Rayleigh's equation when u[z) is analytic 
in some neighbourhood of the point z = Zc in the complex 2:— plane and 
u'{zc) 7^ where u{zc) = c. Their results agree . They show that one 
of the solutions is analytic near Zc and vanishes at it while the other has a 
logarithmic branch point there . 

The results proved here will be described now . In § 2 and § 3 we mainly 
examine the stability of flows for which u'{z) , A{z) and /3{z) are nonneg- 
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ative functions while theorems [3.15] and [3.22] deal with other classes of 
flows . 

Theorem [2.7] is the generalization of the criteria of Rayleigh and Fj0rtoft . 
The only difference from the homogeneous case as above is that the function 
u{z) is replaced by u+{z) , where 2 u^{z) = [ ± ■^^/~K(z) ] . It 

shows that the condition ^ A{z) > for some z ' of Howard and Miles 
is not sufficient for instability . It is also shown that for unstable modes the 
growth rate A; Cj — ^ as A; — ^ oo , and that if /3{z) > everywhere 
then k is bounded above . These prove the predictions of Howard ( 1961 ) . 

Section (3) relates the ffow-curvature T{z) = ^+(2:) / 2u'_{z) to stability. 
Theorem [3.2] gives a pair of bounds k coth [k {z — Zi)] , i = 1,2, which 
T{z) must cross if {k, c) is an unstable mode. In particular T{z) 
must cross the bounds 1 / {z — Zi) if the ffow is unstable . It is shown in 
corollary [3 .4] that if T{z) is bounded above or bounded below then k is 
bounded . Corollary [3 .5] shows that for a fixed upper or lower bound on T , 
ffows with sufficiently small depth {z2 — zi) are stable . 

Theorem [3.6] sets an order in which the bounds must be crossed by T{z) 
in an unstable ffow {u,j3) when u"{z) and j3{z) have no zeros in 
common . If {k , c) is an unstable mode then for some ti < t2 , 
T{ti) = k coth[k {z — Zi)] for each i . Theorem [3 .7] includes the case 
when u"{z) and I3{z) have common zeros . 

Corollary [3.8] shows that if for some Zg in [ zi , 2:2 ] , T is bounded above 
when z < Zg and is bounded below when z > Zg , then k is bounded for 
unstable modes {k , c) . Theorem [3 .13] leads to the same conclusion when 
T is bounded below in z < Zg , bounded above in z > Zg , u'l{z) and 
j3{z) have no zeros in common and u'{zg) 7^ 0. Theorem [3 .14] reffnes 
the result when u'J^{z) and (3{z) have common zeros . 

In § 3.4 stability of ffows is examined when u\z) and A(2;) change 
sign in [ 2;i , ^2 ] in a prescribed way ( see condition B ) . Theorem [ 3 .15 ] shows 
that weaker but nontrivial necessary conditions of instability continue to persist 
showing thereby that some of these ffows are stable . For ffows not satisfying 
condition (B) , no conditions seem to be necessary for instability , suggesting 
that all these ffows are unstable . 

In §3.5 a sufficient condition for the propagation of neutral waves ( Cj = 
0) is obtained. Theorems [3 .17] and [3.18] give necessary conditions for 
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existence of marginally stable modes (umin < c < Wmax) and internal 
gravity waves ( c < Wmin or c > Wmax ) respectively. It is seen that the 
possibility of propagation of marginally stable neutral waves of arbitrarily large 
wave numbers remains open even if the flow-curvature be bounded but the set of 
wave numbers of the internal gravity waves is bounded if either of the functions 
u1.{z) / 2 u^{z) is bounded. 

Upper bounds are obtained in § 3.6 in two different situations , on the phase 
speed Q , for unstable flows not satisfying a Rayleigh-Fj0rtoft type condition . 
Theorem [3 .19] shows that for a flow satisfying condition (A) either U-[z) 
has a point of inflection at Zs with u'J_{z) {u{z) — u{zs)} > for some z 
or Ci — >■ as l^max — ^ where ^/^(z) = [1 - Ag f3{z) / u''^{z)] . Theorem 
[3.22] shows that for a flow with /3{z) < somewhere , either u+{z) has 
a point of inflection at Zg with ^"(z) {u{z) — u{zs) } < for some z or 

Q > as l^max ~^ 1 • 

Propositions [ 3 .16 ] and [ 3 .20 ] together show that the set of complex wave 
velocities of normal modes of an arbitrary flow is bounded . 

In §4 TGE is converted into a quadratic recursion relationship on a 
sequence of smooth functions on [zi, Z2] and solved . For a class of flows 
Theorem [4.1] gives a pair of linearly independent smooth solutions of 
TGE in some neighbourhood of the critical layer z = Zc where u{zc) = Cr , 
when Cj > is sufliciently small . Further if = corollary [4 .3] gets back 
as a limiting case , the smooth solution of Rayleigh's equation that vanishes 
at z = Zc when q = , in agreement with results of Heisenberg and 
Tollmien ( see Drazin & Reid 1981 ) . 

In §4.1 suflicient conditions are suggested for the existence of an unstable 
mode solution to the Taylor-Goldstein boundary value problem . 

2 The basic method 

In this section stability of monotonic flows with nonnegative discriminant will 
be discussed . The basic method is developed in § 2 . 1 . It is then used to obtain 
generalizations of the instability criteria of Rayleigh and Fj0rtoft in §2.2. 
Sharper use of this method will be made in § 3 to obtain deeper results . 



7 



§2.1 Canonical transformations of TGE 

The information on stability obtained from transformations of TGE is stated 
in lemma [2.1]. It is then restricted to a canonical class of transformations to 
obtain propositions [2 .3] and [2 .6] . These have been developed in sufficient 
generality in view of their many applications later in the text . 

Lemma 2 . 1 Let {u, (3) be an unstable flow . Suppose f is a complex- 
valued C^— function. Let f" be integrable and f{z) for every z 
in [zi , Z2] then 

r p\F'\uz = r{ff" + Af)\F\uz. 

J Z\ J Zi 

Further if /m{/^(^)} < for every z in [^1,^2] ; where Im denotes 
the imaginary part then 

Im[ f{z) f"{z) + A{z) f^{z) ] < for some z in (^1,2:2). 

Proof: Substituting F{z) = w{z)/f{z) in (1.1) gives 

[ f\z) F'[z) ] ' + [ f[z) l"(z) + A(z) f{z) ] F[z) = (2.1) 

valid at all points z where f"{z) exists. On multiplying this with F*[z) 
( the complex conjugate of F[z) ) and integrating one gets 

rz2 fZ2 rz2 

/ f\F'\'dz= / (ff" + Af)\F\'dz+ / (fF'F-ydz. 

J Z\ J Zi J Z\ 

Now f'^F'F* = [w'f — wf']{wf~^y is continuous and its derivative is 
integrable in [zi, Z2] . Further f^F'F* vanishes at the boundary points 

zi and Z2 by ( 1 .2 ) . Thus the last integral in the equation above vanishes 
proving the first statement . The lemma now follows from the imaginary part of 
this equation. 

Let the flow-discriminant ^{z) = [ u'^{z) — 4g /3{z) ] . In this and the 
next section we shall mainly focus upon the class of flows {u, (3) satisfying : 

Condition (^4) : (a) u'{z) , P{z) > for every z in [zi, Z2] . 

(6) A(z) has a smooth extensions to some neighbourhood of [^1,^2] in 
R satisfying A (2) > everywhere. 
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Remark 2 .2 It is easily seen that if {u , f3 , k , Cr , Ci ; w) is a solution of the 
equations (1-1) <fe (1-2) then so is { — u , (3 , k , —Cr , Ci ; w* ). Thus either 
both the flows {u ,/3) and { — u,/3) are stable or both are unstable. Therefore 
any result which does not change on u being replaced by — u if proved 
under condition {A) , would remain valid if the hypothesis : u'{z) > 
in (a) above is replaced by u'{z) is monotonic in [zi, Z2] . All the 
results proved under condition (A) here are of this type . 

Suppose the unstable flow {u ,(3) satisfies condition (A) . Let 

1 

u±iz) = uizi) + - [u\z) ± y/X(zj]dz. 



Part (b) of condition (^4) ensures that \J A(z) is smooth so that u±{z) 
are smooth functions and lu'jjyz) = [ w'(^) ± A(z) ] > for every 

z in [zi, ^2] • Further if = then = u{z) . 

Let {k , c) be an unstable mode of the flow {u ,P) . The function 

{u{z) — c } then takes its values in the 3rd and 4th quadrants. 

Let u{z) — c = 1^(2;)— c I exp{iO(z)} (2-2) 

where - tt < e{z) < . (2.3) 

6{z) is then a well defined smooth function of z and a smooth branch of 
\og{u{z) — c} on [zi , Z2] is defined by 

log{w(2;) — c} = log 1^(2;) — c I + i 6{z) . 

On differentiating this equation and taking the imaginary part , one gets 

«'W = (2.4) 

\u{z) — c p 

> for every z in [ 2;i , 2:2 ] • 

Let h{z) be a real-valued , piecewise— function on [zi,Z2] . Let 
Zo be in [ 2:1 , 2:2 ] and let 



f{z) = exp 




(2.5) 
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Clearly f{z) is a C^— function and /" is integrable on [zi,Z2] . Also 
f{z) for any z in [zi, Z2] . 



We write 
where 

so that 



fiz) = \fiz)\ exp[z0(^)/2] 



(f){z) = e{Zo) + Q / 

J z, 



dz 



(2.6) 
(2.7) 



2 Q u'_ (z) 



> 



I u{z) — c p 

It follows now from ( 2 .4 ) and (2.7) respectively that 

> ct>\z) - e\z) = ~^'/^,2 <t>{zo) = e{zo) (2.8) 

I u{z) — Cf 

so that \(f){z) — (l){zo)\ < I ^(2;) — ^(^o) I for every z in [^1,^2]- 



Thus (f){z) lies between 6{zo) and 0{z) . Equations (2.3) and (2.6) now 
give in turn 



and 



-7T < (j){z) < 

Im{f\z)} = If (^)l sin{0(^)} 

< for every z in [zi, Z2] 



(2.9) 
(2.10) 



It is easily calculated from (2.5) that 



mriz) + A{z)f{z) = 



fiz) (2.11) 



{u{z)-c} 

where Vhiz) = [u!l.{z) — 2h{z)u'_{z)'\ . Equation (2.1) then yields 



[f\z)F'{z)]'+i^{h' + h'-k')- 



f\z)F{z) = 



{ u{z) - c } 

On taking the imaginary part of ( 2 .11 ) and using ( 2 .2 ) and ( 2 .6 ) one obtains 

Vh sin (0 - e) 



Im[ff" + Af] = 



{h' + h^-k^) sin - 



u — c 



In view of ( 2 .10 ) and lemma [2.1] it is proved that 
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Proposition 2 .3 Let {u ,13) be an unstable flow satisfying the condition (A) . 
Let h be a real-valued , piecewise—C^ function on [zi, Z2] and let 
zi < Zo < Z2 then for some z in {zi,Z2). 

V,(z) sin{^(z) -e{z) } ^ ^ 

I u{z) — c I sm (f){z) L V / J 

Further suppose [ h'{z) + — A;^ ] < everywhere then the inequality 

above holds at some z which also satisfies (f){z) ^ 6{z) and u{z) ^ u{zo) . 

Remark 2 A If f{z) and <^+{z) are defined by replacing u'_{z) by 
u!^{z) in (2 .5) and (2 .7) respectively then an argument similar to the one 
above yields : 

If — IT < (f)+{z) < then proposition [2 .3] continues to hold provided 
(j){z) is replaced by (/)+{z) and lJh{z) by [u'J.{z) — 2h{z) u!^{z)'\ . 

It will be seen in §3.6 that violation of the bounds on 4>+{z) as above 
implies an upper bound on q . 

Definition 2 .5 We say that two numbers a, b G M , are sign-equivalent and 
write a b if a = Xb for some A > i.e. either both a and b are 
positive or both are equal to zero or both are negative. 

Equations (2.3) and (2.9) show that — 7r<0 — ^<7r so that 

sin I (j)(z) — 9{z) } ~ I (j)[z) — 9{z) } for every z in [^1,^2]- 
Further ( 2 .8 ) implies that 

{ (^{z) - e[z) } ^ - { u[z) - u[zo) } if m ^ e[z) . 

Thus if sin { (t){z) - e{z) } 7^ then 

sin { 0(z) - ^(z) } ~ - {u{z)-u{zo)} . (2.12) 
We now prove : 

Proposition 2.6 Suppose the flow {u,/3) is unstable and condition (A) 
holds then ^{z) > for some z . Let h be a real-valued , piecewise—C^ 
function satisfying [ h\z) + h^{z) — k'^'\ < for every z in [^1,^2] 

11 



and let \Jh{z) = [ u"{z) — 2h{z)u'_{z) ] then 

(a) \Jh{z) takes both negative and positive values in [^1,^2] so that 

U/j,(^s) = o for some Zs in (^i, ^2) • 
(6) U/j(2;) {u{z) — u{zs)} < for some z in{zi,Z2) . 

Proof: Let Zo be any point in [^^1,-22]- Proposition [2.3] shows that 
sin(0 — ^) 7^ for some z so that A(z) > somewhere in [2^1, 2^2] by 
( 2 .8 ). Further (2.9) and (2.12) now give 

\Jh{z) [u{z)—u{zo)] < for some z in {zi,Z2). 

Taking Zq = z\ (resp. Zo = Z2 )here, we see that U/i(z) takes both 
negative (resp. positive) values because u{z) is monotonic. This proves 
part (a) . Now (b) follows on substituting Zo= Zs . 



§ 2 .2 Generalization of Rayleigh , and Fj0rtoft criteria 

Parts (1) and (2) of theorem [2.7] generalize the well known instability 
criteria of Rayleigh and Fj0rtoft while parts (3) and (4) prove the 
predictions of Howard ( 1961 ) . Boundedness of l^ will be proved under 
weaker conditions in § 3 . Corollary [2 .9] shows that the condition of Howard 
and Miles is not sufficient for instability . 

Taking h = in proposition [ 2 .6] gives the parts (1) and (2) below. 

Theorem 2.7 Let the unstable flow {u,/3) satisfy condition {A) and 
let {k , c) be one of its unstable modes then 

(1) u'_l{z) takes both negative and positive values so that 

u^{zs) = for some Zs in {zi,Z2) , 

(2) u'^{z) [ u{z) — u{zs) ] < for some z in {zi , Z2) , 

(3) k"^ Ci < - \ \ max so that the growth rate kci ^ as k ^ 00 , 

(4) // (3{z) > everywhere in [^i, 22] , then, k'^ < [ T'{z) + T'^{z) ] 
for some z in (2:1,2:2) , where T{z) = u'J^{z) / 2u'_{z) . 
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Proof: Only the parts (3) and (4) remain to be proved. It is an easy 
consequence of Howard's semicircle-theorem that Wmin < < Wmax , 
( a result due to Synge ) where u min and u max denote respectively the 
minimum and the maximum of u{z) for z in [zi, Z2] . Let Zq be a 
point satisfying u{zo) = Cr , then ( 2 .2 ) shows that 

exp { i 6{zo) } = — i . 

Thus 6{zo) = —it/ 2 and (j){z) lies between —it/ 2 and 6{z) by (2.8). 
Let sin x = \ sin 9 \ and sin y = \ sin (f) \ where 0<x<y<7T/2. 
Clearly then sin (y — x) = | sin {(j) — 9}\ because both 6{z) and <^{z) are 
either in ( — tt , — tt / 2 ] or in [ — tt /2 , ) . Now 

max [ sin (x) sin (y — x) ~\ = sin^ (?/ / 2) = - sin y tan {y / 2) < - sin y , 

SO that 2 I sin ^ I I sin (0 — ^) I < | sin | for every z in [zi , Z2] . 
Proposition [2 .3] now shows that 

2k^ I u{z) — c I I sin ^ (2;) I < | u^{z) \ for some z and so 

2/c^Cj < \ u'_l I max = max I u'_l{z) I . 

Z€ [zi,Z2] 

This proves part (3) . Part (4) follows from proposition [2 .3] , on substituting 
T{z) for h{z) because Ut = . This completes the proof . 

Remark 2 .8 It is easily seen from the theorem that the criteria of Rayleigh 
and Fj0rtoft remain vaUd , as stated by them , when the fiow-discriminant A(z) 
is a constant or the Richardson number R{z) = g (5 / u''^ is a constant. 

The theorem shows that the flow {2 g z'^ , g z^) is stable though A{z) is 
positive everywhere . 

Corollary 2.9 The condition A (2) > for some z in [zi, Z2] is 
not sufficient to ensure instability. 

^ . . „/ X „/ X u'(z)u"(z) -2g(3'iz) 

Further we have 2u'l{z) = u"{z) H — y . This gives : 

V ^(^) 

Corollary 2 .10 Suppose u'{z) , /3{z) , and A{z) > for every z in 
[zi, Z2] , then 

u'[z) u"[z) ~ for some z in \z\^Z2\ ■ 
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3 Geometry of Taylor— Goldstein equation 

Valuable information on the stability of flows satisfying condition (A) is coded 
in the propositions [2 .3] and [2 .6] of the previous section via the arbitrary 
nature of the point Zq and the function h^z) . This information is decoded in 
this section to bring out into the open the intimate relationship between bounds 
on the flow curvature and stability . 

Theorems [3.2] , [3.6] and [3.7] show that T(z) must cross the 
bounds k coth [ ( z — ) ] , i = 1,2, in a definite order if ( /c , c ) is to be 
an unstable mode. On the other hand corollary [3.8] and theorems [3.13] 
and [3.14] give different kinds of bounds on T{z) each of which ensures 
the boundedness of the wave numbers of the unstable modes . 

In §3.4 the stability of nonmonotonic flows with indefinite discriminant 
is examined. Theorem [3.15] establishes that the class of flows satisfying 
condition ( B ) has stable flows in it . 

In §3.5 necessary conditions are obtained for the propagation of neutral 
waves. Theorems [3 .19] and [3.22] in §3.6 give upper bounds on Cj as 
an alternatives to point of inflection type conditions . 

§3.1 Bounds on the flow-curvature u^{z) / 2u'_{z) 

Definition 3.1 Let £^ = { z e [zi , Z2] \ <(z) ^ , or (3{z) 7^ } , 
and let the flow-curvature T : ^ ^ [ —00, 00 ] be defined by 

Tiz) = [<(z)/2ui(z)] . 

Clearly T is a continuous function in view of u'_ {z) > . It will now be 
proved that in any unstable flow T(z) must cross a pair of bounds . 

Let h{z) = A; coth I A; — + e) } , where e > is sufficiently small , then 
h is a C°°— function on [zi, Z2] and 

h'{z) + h'^{z) — = for every z in [21 , 22] • 

Applying now proposition [2 .6] part(a) one gets 

T{z) — k coth{/c (z — zi + e} > o for some z in ^ . 

In the same way taking h{z) = k coth{/c {z — Z2 — e} gives 

T{z) — k coth{/c {z — Z2 — e} < o for some z in ^ . 
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Now in the limit as e — >^ one obtains the following theorem. 

Theorem 3.2 Let the flow satisfy the condition {A) and let {k^c) 

be an unstable mode, then 

limsup T (t) > k coth { k {z — zi) } for some z in and 
lim inf T (t) < k coth { k {z — Z2) } for some z in 3^ . 

t^z 

Remark 3 .3 Suppose u+iz) and (3{z) have no common zeros in 
[z\ ^ Z2] , then T is continuous on [zi , Z2] and the left hand sides of the 
inequalities in the theorem reduce to T (z) . 

Theorem [3.2] has some interesting consequences. Using the inequality 
coth (t) > 1 if t > one obtains : 

Corollary 3.4 Suppose the hypothesis of theorem [3.2] holds. Let T (z) 
be bounded above or bounded below , then 



k < min 



sup{T(^) } , - inf {T(^) } 



< 00 



SO that the set of wave numbers of unstable modes of {u,/3) is bounded . 

In § 3.2 and § 3.3 boundedness of k will be proved under weaker conditions 
Using k coth {kt) > 1 / 1 for t > one obtains : 

Corollary 3.5 Suppose the hypothesis of corollary [3.4] holds, then 
(a) < T (z) for some z in {zi , Z2) and 

(6) — ^ — < — T (z) for some z in {zi , Z2) so that 

Z2 — z 

^ ' sup { T (z) } , - inf { T (z) } < 00 

Z€^ 



(c) < min 

Z2 - Zi 



This shows that given a fixed upper or lower bound on T(2;) fiows with 
sufficiently small depth {z2 — zi) are stable . Thus we see that instability 
needs sufficient room to manifest itself. 
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§ 3 .2 The main instability criterion 

We are now in a position to prove the main result on instability . 

Theorem 3.6 Let {u,/3) be an unstable flow satisfying the condition {A). 
Suppose u!l_{z) and (3 {z) , have no common zeros in [2:1,2:2] then for some 
ti, t2 in [zi , Z2] , ti<t2, 



and 



T(ti) = k coth [k [ti - zi)] > 
T{t2) = k coth [k{t2 - Z2)] < 



1 



ti - Zi 

1 



t2 - Z2 

Proof : Remark [3 .3] shows that T , in this case is a continuous function 
into [ — 00, 00 ] . 

Let 



and 



= min {z \ T{z) > k coth [k {z- Zi)]) 
t2 = max { z I T(2:) < k coth [ A; (2: — 2:2) ] | . 

Zi<Z<Z2 



Theorem [3.2] shows that ti and t2 are well defined , and clearly 

ti ^ t2 and T{ti) = k coth [ k {t^ - z,) ] for i = 1, 2 . 

We prove now that ti < t2 ■ If possible suppose ^2 < • Let Zs be any 
point in (^2,^1) , and Tg = T(2:s) . Let e > be a small number and 
let h : [zi , Z2] — ^ M be defined by 



h{z) = 



k coth [ A; ( 2: — 2:1 + e ) ] if 

h{z) if 

hiz) if 

k coth [ A; ( 2: — 2:2 — e ) ] if 



z e [zi, Zs - e] 

z G [ 2s - e , 2s 

z G [ 2s , 2s + e 

z G [ 2s + e , 22 ] 



where the graphs of liiz) for i = 1,2 are obtained by joining the point 
( 2s , Ts ) with the points ( 2s =F e , ^ coth [ A; ( 2s — 2j ) ] ) respectively , by 
line segments . When e is sufficiently small it is easily checked that , 



1. h is a piecewise — function, and h{zs) = T(2s) ; 

2. /i' + /i^ — A;^ = if 2 is not in [ Zs — e, Zs + e] ; 
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3. k coth [ k {t2 — Z2) ] < h{z) < k coth [ k {ti — Zi) ] for every z in 
[zs-e, Zs + e] ; 

4. for any fixed m < , h'{z) < m for every z in [ — e , + e ] ; 

5. h' -\- h"^ — k"^ < for every 2; in [zi , Z2] ; 

6. T'(2;) — h'{z) > for every 2; in [ Zg — e , Zg -\- e ] ; 

7. T(z) < h{z) li z < Zs-e and T(^) > h{z) if z > + e ; 

8. T(^) = h{z) only if z = Zs and T(^) - h{z) ~ (z - z^) . 

It follows that [ T(z) - ] [ ) ] > for every z in [zi , Z2] . 

This contradicts proposition [2.6]. Thus < ^2 , and the proof is complete . 

With a little more care one can include the case when u'_l{z) and (3{z) 
have common zeros . We state the result and briefly indicate its proof . 

Theorem 3.7 Let {u,/3) be an unstable flow satisfying condition {A), then 
for some ti and t2 in ^ , t\ < t2 , 

(A) limsup T(2;) > k coth \k{t\ — z\) 1 > — ^ — and 

z^ti ti — Zi 

(B) liminf T(^) < k coth [ k {t2 - Z2) 1 < — ^- — 

z^t2 t2 — Z2 

and one of the following conditions (a) or (b) hold: 

(a) ti < t2 ; 

(b) ti = t2 = to { say ) and one of the following conditions hold : 

(1) > k coth [ k (ti — zi) ] where L+ = limsup T{z) . 

z^to-O 

(2) R_ < A; coth [k {t2- Z2) ] where R_ = liminf T{z) . 

(3) L+ > R_ . 

T(2^) — To 

(4) L+ = R_ = To (say) and lim inf — — = — 00 . 

z^to Z — to 

Proof: 

Let ti= ini { z e ^ \ limsupT(t) > k coth [ k {z - zi) ] } 

z 
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and t2 = sup { z G ^ | liminf T(t) < k coth [ k {z — Z2) ] } . 

It is clear that inequahties (A) and (B) above hold . 

Suppose both conditions (a) and (b) do not hold . Clearly t2 < h . If 
t2 = h then k coth [ /c (t2 — ^2) ] < L+ < R_ < A; coth [ A; (ti — 2;i) ] and 
ti is not in ^ . Let Zs = [h -\- 12] / 2 and let 

^ _ f T{zs) if t2 < ti and Zs is in ^ 

^ ~ \ ( L+ + R_ ) / 2 if t2 = h or Zs is not in ^ . 

Let the function h be defined as in the proof of theorem [3.6] above . It is 
not difficult to check , as in the proof of the previous theorem , that for sufficiently 
small e 

[ T{z) — h{z) ] [ u{z) — u{za) ] > for every z 'm.[zi^ Z2\ ■ 

This contradicts proposition [2.6]. Thus condition (a) or (b) above must 
hold. 

For any ^3,^1,^2 in [ ^1 , ^2 ] , < ^2 , either is in [ zi , ] 

or t2 is in [ , ^2 ] • From parts (A) and (B) of the theorem it follows 
easily that : 

Corollary 3 .8 Let {u, (3) be an unstable flow satisfying condition {A) , 
and let Zg be a point in [zi, Z2] . Suppose T{z) < Tq in z < Zg 
and T{z) > — Tq in z > Zs for some > , then k < Tq and 
the set of wave numbers of unstable modes of {u ,/3) is bounded . 

§ 3 .3 Bounds on the wave numbers of unstable modes . 

Corollary [3.8] gives conditions on T(z) so that the unstable modes have 
bounded wave numbers . In this subsection different sets of conditions will be 
described that ensure the same . The following lemma generalizes a deft piece 
of estimation in Drazin & Howard ( 1966 ) . Let Zs ^ [zi^ Z2\ and let 

u{zi) < 2Cr — u{Zs) < U{Z2) j 

2 Cj, — u{zs) < u{zi) ; 

2Cr — u{Zs) > U{Z2) . 



2Cr — u{Zs) 



Uq 



U{Zi) 
U{Z2) 



if 
if 
if 
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Lemma 3.9 Let Zg [zi, Z2] and Uq be as above. Let Zq be a point such 
that u{zo) = Uq. Let 6{z) and (f){z) be as in proposition [2 .3] , then 
vi^) ^ 1 for every z in [zi,Z2] where 

^^^^ _ u{z) - u{zs) sin {0(z) - e{z)} 



u{z) — c I ' sin (f){z) 

Proof: let Zc be a point such that u{zc) = Cr . For any t , we write 
Ut and 6t for u{zt) and 0{zt) . We assume first that Us < Cr , then 

Zs ^ Zc < Zo and 6s < — 7r/2 < Oq < —ir — Os- 

Case(l) : zi < z < Zg , then 

u{z) < Us and < (f){z) - e{z) < tt + e{z) < - Os < it , 
so that Ci cot(0 — 6) > — Ci cot(^) = Us — Cr . 

{u — Us) 



Thus r]{z) 



II r \ COS 9 — sinO cot(0 — 0) } 

{u - Us) 



{U — Cr) — Ci COt(0 — 6) 
^ {Us - U) ^ ^ 

~ {Cr - w) + {Us - Cr) 

Case(2) : Zs < z < Zq , then 

u — Us>0, < (f) — 6 < TT and sin0 <0 so that r]{z) <0 

Case (3) : Zo < z < Z2 , then 

(u-Us) > and > - ^ > ^0 = - - > - tt 

and so q cot(0 — 0) < q cot(— tt — 6s) = {us — Cr) . 

U - Us) 



Thus r]{z) 



(11 — Cr) — Ci cot(0 — 6) 



^ {U - Us) ^ ^ 

~ {u — Cr) — {Ug — Cr) 

We have proved that if Us < Cr then 77(2;) < 1 for every z . A similar 
argument applies when Us > Cr and completes the proof of the lemma . 

Proposition 3 .10 Let {u , (3) be an unstable flow satisfying condition {A) . 
Let h be a piecewise—C^ function on [zi, Z2] ■ Let Zs be in [zi, Z2] 
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and \Jfi{zs) = . Let Kh{z) = — \Jh{z) / { u{z) — u{zs) } whenever 
u{z) 7^ u{zs) and suppose Kh{z) > . Let k'^ be sufficiently large 
so that [ h'{z) + — /c^ ] < for every z in [zi , Z2] , then 

k'^ < sup { h\z) + h^{z) + Kh{z) \ u{z) ^ u{zs) } 

z 

and so the set of wave numbers of the unstable modes is bounded if 
K/i(2;) is bounded on the set { u[z) 7^ u[zs) }. 

Proof: Let Zq be as in lemma [3 .9] . From proposition [2 .3] we have 

¥^h{z) . ri{z) < h' (z) h'^ (z) — k'^ for some z in [^1,^2]- 
The proposition now follows from lemma [3 .9 



Remark 3 .11 Suppose u'{zs) 7^ 0, then K(2;) has a continuous 
extension to [zi, Z2] , and so it is bounded . 

Taking h = in this proposition we obtain a generalization of a result 
due to Drazin &; Howard ( 1966 ) . 

Corollary 3.12 Let Ko(z) = —u'l{z) / { u{z) — u{zs) } > 0, whenever 
u{z) 7^ u{zs) , then k^ < K.q{z) for some z . 

We use this proposition to obtain the following intrinsic criterion for the 
boundedness of the wave numbers of unstable modes . 

Theorem 3.13 Let {u^(3) be an unstable flow satisfying condition {A). 
Suppose u'l{z) and (3{z) have no common zeros. Let zi < Zg < Z2 
and let u'(zs) 7^ 0. Let T (2;) > Ti if z < Zs and 
T{z) < T2 if z > Zg for some constants Ti and T2 , then 
the set of wave numbers of the unstable modes is bounded . 

Proof: It is clear from the hypothesis that T(zs) is finite. Let 
Tg = T(zs) and let e > be a small number. Let h be the piecewise 
linear function on [^1, ^^2] ) whose graph is obtained by joining the points 
( , Ti ) , ( - e , Ti ) , ( , Ts ) , ( + e , T2 ) and ( ^2 , T2 ) . 

Clearly h{zs) = T(zs) . When e is sufficiently small and u{z) 7^ u{zs) , 
we have 

^ > so that K(z) = —-- — > • 

u{z)-u{Zs) u{z)-u{Zs) 
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Further by the remark above , K(2;) is bounded above . The theorem now 
follows from proposition [3 .10] . 

With some more effort one can include the case when ^+(2;) and /3{z) 
have common zeros . We state the result . 

Theorem 3.14 Let {u.(3)he an unstable flow satisfying condition {A) . 
Let zi < Zs < Z2 . Let u'{zs) ^ . Let T (z) satisfy 

T (z) > Ti if z<Zs and T{z) < T2 if z> Zs 
for some constants Ti and T2 . 

(b) L_ > R+ , where L_ = liminf T(z) and R+ = limsup T(z) 



Then the set of wave numbers of the unstable modes is bounded . 

§ 3 .4 Flows with indefinite velocity and discriminant . 

In this subsection the stability of a class of flows , in which u\z) and A(z) 
change sign in [zi, Z2] in a prescribed way will be discussed . It is shown 
that nontrivial necessary conditions of instability continue to persist , though in 
a weaker form , showing thereby that some of these flows are stable . This class 
of flows is described by the following condition . 

Condition(B): Allow ( u ,/?) satisfies this condition if for some Si , S2 in 

[ Zi , 2:2 ] , Sl < S2 , 

(1) P{z) > for every z in [zi, Z2] . 

(2) A(z) has an extension to some neighbourhood of [zi, 22] satisfying 
A(z) < if z is in [§1,52] and A(z) > otherwise. 

(3) u{z) is monotonic in each of the intervals, [zi. si] and [s2, ^2] • 
Thus A{si) = , for i = 1,2 ; u'{z) is unrestricted in [ si,S2 ] , 

and its signs in Si = [zi , si] and ^2 = [s2, Z2] could be different . For 
i = 1 , 2 let 



Z—^Zs+O 



(c) If L_ = = Tg (say) , then limsup 



T{z) - T, 



7^ 00 . 



Z- Zs 




1 
1 



if u'{z) > for every z in Si ; 
if u'{z) < for every z in Si . 
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1 

Let u±{z) = u{zi) + ^ / [ u'{z) ± p{z)]dz , 

where r)(z) - i V'^(^) ^(^) > and z is in 5, ; 

wiftere P[z) - <^ ^ ^^^^ ^ ^ 

This extends the definition of u±iz) given earher after remark [2.2]. 
u±{z) are C^-functions in [^^1,^2] , and have continuous second derivatives 
except possibly at the points Si and S2 ■ We now prove : 

Theorem 3 .15 Let {u,/3) be an unstable flow satisfying condition (B) , then 
(a) {—Y (7iu'l{z) < for some i = 1 or 2 and some z in Si . 

(6) // A'{si - 0) ^ A'{s2 + 0) then u+{z) has a point of 
inflection at some point Zg outside [si, S2], satisfying A{zs) < 0. 

(c) // Si = S2 5 ^G-t one of the conditions (1) or (2) hold : 

(1) L+ < R_ where L+ = limsup T{z) and R_ = Hminf T(z) , 

T ( T 

(2) L+ = R_ = To (say) and Hminf 7^ —00, 

z^ta Z — to 

then one of the conditions (3) or (4) is satisfied : 

(3) limsup T (t) > k coih. {k[z — Zl)} for some z in [2:1,81], 

t^z 

(4) liminf T (i) < /c coth { A; ( 2; — 2:2 )} /or some z in [82,2:2], 

t^z 

so that if (— )*"^^T(^) is bounded above in Si for each i, then 



k < max 



sup{T(z) } , - inf {T(z) } 



< 00 



Proof: Let Zq = [si + S2]/2. Let (/){z) be as in (2.7), with 

u±{z) as defined above. Clearly (t>'iz) = 0'{z) in [si,S2]. It 

follows that 0(2:) = d{z) in [si, S2]. Further if z is in Si and 
(f)[z) ^ Biz), then [ (f){z) - e{z) ] ~ Oi [ see definition 2 .5 ] , 

because (t>'{z) — 9'{z) ~ 2u'_{z) — u'{z) ~ —<Ji. 

It follows from proposition [2.3] and equation (2.9), that for any 

piecewise-C^ function h{z) , satisfying [ h'{z) + /i^(2:) — /c^ ] < and 
for some i = 1 or 2 

(-)' (7i Uft(2:) < for some z m. Si. (3.1) 
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Tn = 



Part (a) follows from this on taking h = . 

Under the conditions of part (b) , ai u^{z) — ^ — oo as 2; — ^ si — , 
and a2u'l{z) — ^ 00 as 2; — ^ S2 + . Now part (b) follows from (a) . To 
prove part (c) we may assume that if si = S2 then 
k coth [ /c (§2 — -2^2) ] < L+ < R_ < /c coth [ A; (si — Zl) ] . Let 

T(zo) if Si < S2 and Zq is in ^ 

[ L+ + R_ ] / 2 if Si = S2 or Zo is not in 3^ . 
Let e > be a small number . Let h{z) be the function defined in 
the proof of theorem [3.6] with Zq and Tq replacing Zg and Tg 
respectively . If possible suppose neither of the conditions (3) and (4) hold . 
It is clear then that for sufficiently small e , 

T{z) - h{z) < in and T{z) - h{z) > in S2 . 

Further aiu'{z) > for every z outside [si, 52]- It follows that 
(— cr^ [ u1{z) — 2h(z)u'_{z) ] > for each i = 1 and 2 and every z 
in Si . This contradicts (3.1) and proves part (c) . 

§ 3 .5 Propagation of neutral modes 

A normal mode ( /c , c ) of a flow {u ^(3) is called a neutral mode if q = . 
It is called a marginally stable mode if c is in [u mm , u max ] • These modes 
are important . In the space of normal modes of flows where a point looks like 
{u , P , k , c) satisfying (1.1) and (1.2), it is the marginally stable modes 
that separate stability from instability . 

A neutral mode with c outside [ w min , u max ] is called an 
internal gravity wave. Lemma [2.1] with f(z) = {u{z) — c} then gives 

- cf\F'\''dz = r [-k'^u - cY + g(3{z)] \F{z)\''dz . 

J Zl J Zl 

This shows as observed in Drazin & Howard ( 1966 ) that no internal gravity 
waves can exist if (3{z) is negative everywhere . Otherwise using the well 
known inequality that if a C^— function ip{z) satisfies ip{zi) = = ip{z2) 
then 

/ \Hz)\^dz < ^ ' / \,lj\z)\^dz , (3.2) 

Jzi Jzi 

one obtains \u{z) — c\'^^ < g(3inax{z2 — zi)"^ / n'^ . It follows that : 
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Proposition 3 .16 Let {k,c) be a neutral mode of a flow {u,f3) 
and suppose /3{z) > for some z , then 

Wmin - a < C < Wmax + Ci where a = VS'^max [Z2 - Zi] / 7T . 



Let c G M. Let P{zs) = /3'{zs) = = u"{zs) whenever u{zs) = c. 

T^/ \ u"(z) g (3iz) ^ . I 1 . r 1 

Let Kfz) = p — ^ — I — p — T-z- , be Integra Die m zi , Z2\ ■ 

{u[z) — c\ {u[z) — 

This is the case if c is outside [wmin, Wmax] otherwise the condition 
u'{zs) 7^ whenever u{zs) = c is sufficient to ensure this . Suppose K(z) > 
/ { — ^1 }^ whenever u{z) ^ c. Equations (1.1) and (1.2) then 
constitute a regular Sturm-Liouville problem , and so the flow {u , (5) admits a 
neutral mode ( fc^ , c ) , where 



k1 = — min 



I -K(t)'')dz I j\''dz^ >0 



the minimum being taken over all functions such that and (J)' are 
square integrable and 0(^i) = = 0(2:2) . Further if (A;, c) is any neutral 
mode then < . We now prove : 

Theorem 3.17 Let {k,c) be a smooth neutral mode for {u , P) satisfying 
condition (A) . For some Zs ^ [ ^1 , ^2 ] let u (zs) = c , (3{zs) = = 
u1{zs) and u'{zs) 7^ 0, then one of the following conditions hold: 

(1) limsup T(t) > k coth [ /c {z — zi) ] > — ^ — , for some z < Zg . 

t^z Z Z\ 

(2) liminf T(t) < k coth [ A; {z — Z2) ] < , for some z > Zs . 

(3) Condition (6) of theorem [S .7 ] holds with ti = t2 = to = Zg . 
Proof: u'_{zs) = because (3{zs) = 0. Let h{z) be piecewise— C-^ . 



Let f{z) = exp 




{u{z)-c} 



+ h{z) > dz 



f{z) is then smooth and positive everywhere . Equation ( 2 . 1 ) then shows 
that for every z in [-21,-^2], 



[f\z)F'{z)]' + 



F{z) = 
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Multiplication by F*{z) and integration as in lemma [2.1] yields 



{u{z) - u{Zs) } 



< { h'{z) -\- h {z) — k } for some z & {zi , Z2) 



The theorem follows from this by an argument similar to that used to prove 
theorem [3.7] . We next prove : 

Theorem 3.18 Let {u,P) satisfy condition (A) . Let {k,c) be a 
smooth neutral mode and let c be outside [wmin ? i^max] ? then 



u'![z, 



max 



(1) — , , > k coth \ k iz — zi) ^ for some z if c > u 

(2) 2'^' [ ) ~ ^ ^^^^ [ k {z — Z2) ] for some z if c < Wmin- 

(3) Parts (1) and (2) hold if u+{z) and U-{z) are interchanged. 

(4) // either of [ u![{z) / 2 u^{z) ] is bounded then k is bounded. 



Proof: Let f{z) = exp 



+ h{z) > dz 



in turn 



{u{z)-c} 

Parts (1) and (2) follow on taking the — ve sign above by an argument 
similar to that used to prove theorem [3 .6] while part (3) follows on taking the 
+ ve sign . Part ( 4 ) follows from ( 1 ) , ( 2 ) and ( 3 ) . 



§ 3 .6 Bounds on the phase speed Cj of unstable modes . 

This section deals with two situations where upper bounds on the phase speed 
Ci are obtained as alternatives to Rayleigh-Fj0rtoft type conditions . For flows 

satisfying condition (A), theorem [3 .19] is obtained if U-{z) is replaced by 
u+{z) in (2.7). Theorem [ 3 .22 ] is about flows in which /3{z) is negative 

somewhere. In both cases an upper bound on q is obtained if (2.9) does 

not hold . 

Let A = TT + = tan~-^ [ Q / (c^ — Umin) ] , 

and B = - 9{z2) = tan~^ [ q / (wmax-c^) ] . 
Clearly < A , B < ir / 2 . Let r = |[wmax — Wmin] then Q < r . 
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It follows from (2.4) that 

e{z,) - e{z,) = n-A-B = c, ru\z)dz_ 

Jzr I W - C I ^ 

Further min{A,B} > tan"^ f — )> — tan"^ (1/2) (3.4) 

V 2 r / r 

because tan~^ x > [x / Xq) tan"-^ if 0<x<Xo<7r/2. 



and A + B = tan"^ 



2rci 



> tan-^ {2ci/r) > — tan"^ (2). (3.5) 

r 

Theorem 3.19 Let {k ^ c) he an unstable mode for the flow {u,(3) 
satisfying condition {A) . Let y{z) = A(2;) / u'{z) and = max i/{z) , 

z 

then one of the following conditions (a) or (b) holds 

(a) u'J_{z) takes both negative and positive values so that u"{zs) = 
for some Zs in {zi, Z2) , and u"{z) [ u{z) — u{zs) ] > for 
some z in (zi , Z2) . 

(6) c^ < —- — iJ^^"^^^ — 1 — _i /o^ 1 ^^"^ Q ^ as r z/^ ^ . 
[ tan ^( 1 / 2 ) + Vrn tan ( 2 ) J 



Proof 



2u' (z) 

Let (t>+(z) = Oizo) + Ci —. — -— ^ 77^ dz and suppose the 

Jzo I U{Z) - C |2 

condition: — tt < 4)+{z) < for every Zo,z in [zi, Z2] , holds. 
In this case it is easily seen that sm{ (j)+{z) — 0{z)} ~ {u{z) — u{zo)}. 
Remark [2.4] then leads to the result in proposition [ 2 .3] with u+{z) and 

U-{z) interchanged and replaced by 0+ . The arguments used to prove 
parts (1) and (2) of theorem [2 .7] now prove the part (a) above. 

On the other hand suppose for some Zo, z in [^1,^2], the condition : 

— TT < ^+(^) <0 does not hold, then 

c^ r ^""^/"^ > min(^, 5) + eiz2) - e{z{) 

Jzi \U - C\^ 

SO that Umi'^ — A — B) > min {A, B) 

because of (3.3). Part (b) now follows from (3.4) and (3.5). This 
completes the proof of the theorem . 
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The flows {u ,/3) with /3{z) < somewhere wifl be considered now . In 
this case z/^ > 1 . All these flows are believed to be unstable . We now prove 
that for such a flow the complex wave velocities of unstable modes is bounded . 

Proposition 3.20 Let {k,c) be an unstable mode of a flow {u,P) and 
P{z) < somewhere. Let = -\- [ — fif/^lmax (^2 — ^i)^/7r^5 then 

(1) ^min < Cr < Wmax and 

(2) [ Cr - K^max + ^min ) ] ^ + cf < . 

Proof: Taking /(z) = {u{z) — c} in lemma [2.1] and proceeding as in 
the proof of Howard's semicircle theorem ( see Howard 1961 ), one obtains 

(1) / \u{z) - Cr}{\F'\'^ + k^\F\'^)dz = and 



(2) { [ - i( nn,ax + Wmin ) ] ' + cf - T^} H { \ F ' \^ ^ F \^ ) d . 

< - r gP\F\^dz < [-gPUs. r\F\Uz. 



' Zl J Zl 

The proposition now follows from these and (3.2) . 

Remark 3.21 It is clear from propositions [3.16] and [3.20] together 
with Howard's semicircle theorem that for an arbitrary flow {u ,/3) , the set 
of complex wave velocities of its normal modes is bounded . 

Theorem 3 .22 Let {k,c) be an unstable mode for the flow 
{u,l3) . Let u' and A be nonnegative functions and let (3{z) 
be negative somewhere. Let = v{z^ where u{zc) = Cr and let 
a = max { i/'{z) / u'{z) } then one of the conditions (a) or (b) holds. 



(a) Parts (1) and (2) of theorem [2 .7 ] hold 

{Urn - I) 7^ To 

[ tan-i(ro/2r) + ( z^^ - 1 ) tan-i ( 2r J r ) ] 



(6) Ci < — T- ; r r 7- ; ttUd 



[ tan-i(ro/2r) + ( z^^ - 1 ) tan'^ ( 2r J r ) ] ' 
so that Ci ^ as a To ( z^m — 1 ) — ^ . 

Proof: Let (f)(z) be as in (2.7). Suppose the condition: — tt < (J){z) <0 
for every Zq ^ z in [^1,^2], holds . In this case the same argument that is 
used to prove parts (1) and (2) of theorem [2 .7] , proves part (a) . 
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On the other hand suppose for some Zq and z in [zi , Z2] , the 
condition : — tt < (f){z) < does not hold 

then Q / ^ ^- < - mm{A, B) (3.6) 

Jz, \u - c\^ 

so that (z/^-l)(7r-^-B) > min {A, B) . 
The definition of the angles A and B given above , gives in this case 

mm{A,B} > tan"^ f^)> tan"^ (ro/2r) 

V 2 r / To 

and A -\- B > tan"^ ( 2 q / r ) > — tan"^ ( 2 To / r ) , 

SO that ( — 1 ) TJ" > — [ tan~-^( To / 2r) + ivrn — ^) tan"-^ i2ro / r) ] . 

The first inequality in part (b) follows now from this equation . To obtain 
the second inequality, sharper estimation of the integral in (3.6) is needed. 
It may be assumed that v^. < 1 . Let s = (1 — v^) j a then 

v{z) — l<z/c — l + as = if \u{z) — c\ < s 
1 c 

Let E = -7T — tan~^ ( ~ ) • Equation (3.6) then gives 

{um - l)i7T - A - B - 2E) > mm{A,B) 

so that ( z/^ - l){n-2E) > min ( ^ , 5 ) + ( - 1){A + B) . 
Further n — 2E = 2tan~^(ci/s) < 2ci/s. It follows that 

(z/^ - l)a > ""'^ [ tan-i(r«/2r) + (z/^, - 1 ) tan^^ ( 2 r, / r ) ] . 

2 To 

The second inequality in part (b) of the theorem follows from this . 

Remark 3 .23 It is clear from the proofs that in both the theorems above, 
results similar to those of § 3 will hold if the condition in part (b) does not hold . 

4 Fundamental solutions of T G E 

In this section TGE will be solved for a class of flows . It is seen here in 
a partially discreet form as a quadratic recursion relationship on a sequence 
of smooth functions . A pair of linearly independent solutions is obtained in 
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some neighbourhood of the critical layer z = Zc where u{z^ = c-r i in 
theorem [4.1] when q > is sufficiently small . To prove the smoothness 
of the solutions it is assumed that [ g /3{z) / u''^{z) ] is a constant. It is 
very likely that the results hold without this . Corollary [4 .3] gives a smooth 
solution of Rayleigh's equation which vanishes at z = Zc when q = . 

Let U be an open subinterval of {zi, Z2) . Let Zc be a point in U . 
Let u{z) , /3{z) , and v{z) be C°°— functions on U . For z in U let 

w{z) = exp J^^ v{z) dz satisfy (1.1), then 



v'{z) + v\z) - 



+ 



= 



{u{z) — c} ' {u{z) — c}^ 
For n > let an{z) be C°°— functions on U and let 



{u{z) - c} 



+ an{z)Y''{z), where y(^) = 



n=0 



Suppose these series converge uniformly on some open subset J of U 
Substituting these values of v{z) in turn , in the equation above , one gets 



— u^{z) + 2u!j_{z)ao 
{u{z) - c} 

2w4(z) 



+ 



1 + 



00 



E 

n=l 



n+ 1 + 



u'{z) 



I CLn+1 + CLn 



2u!^{z) 
u'{z) 



I ai + -\-al - k' 



u"{z) ^ 

^ On TT^ + 2^ 

j=0 



u'{z) 



+ 



This proves part (a) of the following : 

Theorem 4.1 Let iu,(3) be as above. Let (3{z) , A{z) > and 
u'{z) > for every z in U . Let {k,c) be given and u{zc) = Cr ■ 



Let 



wi{z) = exp 



and 

where ao{z) = 
an+i{z) = 



W2{z) = exp 




n=0 

00 



2u;(^) 

1 



ai(2;) = 



aQ[z} 



ahz 



)} 



1 + 2{u:,iz)/u'iz)} 




(n + 1) + 2{</w'} 









r «" 1 


— ajj + na„ < 







j=0 
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for n > 1 and bn{z) is obtained from an{z) on interchanging 
u+{z) and U-{z) . 

(a) Suppose the series ^ an{z)Y^{z) [resp. ^ bn{z) (z) ] converges 
uniformly in some neighbourhood J of Zc, then wi{z) [resp. W2{z) ] 
is a smooth solution of TGE in J . 

(6) Let Richardson number { g f3{z) / u''^{z) } be a constant in U . Let 
S > 1 + A;^ be a constant . Suppose for every z in U , 
r u"{z) \ 



1 



< S-'"'""^ for every j > and 



dz^ { u'{z) 
(2) IfloC-s;)! ( resp. \bQ{z) \ ) < S, then 

for sufficiently small q > , ^ CLni^z) Y^{z) [ resp. ^bn{z) Y^{z) ] 
converges absolutely and uniformly in some neighbourhood J of Zc- 

Proof: Under the hypothesis of part (b) , u!j_{z) are constant multiples of 
u'{z). It is then clear that a„ and 6„ are polynomials in [u"{z) / u'{z)] 

and its derivatives with constant coefficients . 

d^ ( u"{z) 

dz^ \ u'{z) 

ring of polynomial functions of ? ? • • • with real coefficients . For a multi- 
index m = ( mo , mi , . . . ) with only finite number of nonzero entries , let 
deg (m) = Y^^q U + 1) and t"^ = U^o C • 

For q{z) = ^ in , let deg {q) = max { deg (m) | 7^ } 

and \\q\\s = I I S ^"^^ . It is easily checked that if q, Qi , q2 ^ R 
and A G M , then 



Let tj{z) = -J— \ —777^ } for j > . Let i? = R [to , , • • • ] be the 



(1) sup\q{z)\ < II g lis ; (2) ||gi + g2||s < lki ||s + 11^2 ||s ; 

z 

(3) II gig2 lis < Iki lis Ik2 lis ; (4) || Ag||s = | A| ||g||s and 
(5) II g' lis < s deg{(i) ||g||s . 

The last part follows because on differentiating a monomial of degree n one 
obtains a sum of at most n monomials of degree ( n + 1 ) . 

It is clear from the hypothesis that fln , &n ^ Li for every n > and 
deg {an) = deg {bn) = ( n + 1 ). Further ||ao||s ^ S and ||ai||s ^ 3s^. 
We now prove that Ha? ||s ^ for every j > 0. Inductively we 

assume that this result holds for j < n ^ then 
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Q^n s ^ 



1 

n 



a 



n-l 



+ 



^ n— 1 

n ^ 

j=0 



,n+l 



< 3" s 



n Q n+1 



SO that 



sup I an{z) 



< 3s 



It follows that 
in U , and 



lim sup 

X^^o '^ni^) converges absolutely and uniformly for z 
I y I < r , for any r < 1 / [ 3 S ] . 
A similar argument with replaced by 6„ gives the same result about 

the convergence of Yl^=o bn{z) V^^z) . Part (b) now follows from the fact 
that for sufficiently small q > , the set { | ^(2;) | < r } is a closed 
neighbourhood of Zc in [zi , Z2] . 



Remark 4.2 (1) The function u{z) together with its two derivatives can 
be uniformly approximated by a function q{z) = J^^fexp { J^^^ p [z] dz} dz] 
where p {z) is a polynomial in z . Clearly then [ q"{z) / q'{z) ] together 
with all its derivatives is uniformly bounded . 

(2) Suppose Ci = = (3{zc) = u'l{zc) and T(z) has a smooth extension 
to [zi,Z2],then W2{z) is a smooth solution of TGE . 

When (3 = 0, we have u+{z) = u{z) , and U-{z) = 0. This leads 
to the regular solution of Rayleigh's equation. 



Corollary 4 .3 Suppose f3{z) = and u'{z) > for every 
U. Let ai{z) = / 3 , and for n > 1 



1 



— a' + nttr 



Let 



n + 3 

w{z) = {u{z) — c} exp 



n-l 



^ n—j 



Suppose for some S > 1 + /c^ / 3 , 



J > and every z 



' u"{z) \ 



d^ 

~dz^ 



u'iz 



< S^^^ for every 



in U . Let Q > be sufficiently small then 
the series Yl'^=o '^n{z) ^^{z) converges absolutely and uniformly in some 
neighbourhood J of Zc where u{zc) = Cr and w{z) is a smooth 
solution of Rayleigh's equation in J . 
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§4.1 Existence of unstable flows. 

So many parallel flows are believed to be unstable because they are not seen to 
persist , yet very few can be proved to be so . The only known example of an 
unstable parallel flow with rigid and finite boundaries is a piecewise linear fiow 
called a shear- layer ( see Drazin & Reid 1981 ) . The difficulty is that no general 
result for the existence of solutions of Taylor-Goldstein boundary value problem 
is known . A good beginning in this direction has been made ( see Friedlander 
& Howard 1998 ; Friedlander 2001) but much remains to be done here. An 
existence theorem will provide deep insight into the subject . We expect : 

(1) Let {u,P) be a flow satisfying condition (A) and let T(z) be 
bounded . Suppose for every k'^ < | T | max , every c inside Howard's 
semicircle and every piecewise— function f{z) satisfying /m[/^] < 
everywhere , the condition : Im [ //" -|- Af'^ ] < ^^r some z , holds then 
the flow is unstable . 

(2) Suppose for some to < h < t2 in [zi, Z2] , A(to) < 0, A{t2) < 
and A(/;i) > then the flow is unstable . 

(3) If P{z) < for some z then the flow is unstable . 
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